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In this paper, we present a method for trajectory generation and adaptive full-state 
feedback control to facilitate spacecraft formation flying near the Sun-Earth L2 Lagrange 
point. Specifically, the dynamics of a spacecraft in the neighborhood of a Halo orbit 
reveals that there exist quasi-periodic orbits surrounding the Halo orbit. Thus, a spacecraft 
formation is created by placing a leader spacecraft on a desired Halo orbit and placing 
follower spacecraft on desired quasi-periodic orbits. To produce a formation maintenance 
controller, we first develop the nonlinear dynamics of a follower spacecraft relative to the 
leader spacecraft. We assume that the leader spacecraft is on a desired Halo orbit trajectory 
and the follower spacecraft is to track a desired quasi-periodic orbit surrounding the Halo 
orbit. Then, we design an adaptive, full-state feedback position tracking controller for 
the follower spacecraft providing an adaptive compensation for the unknown mass of the 
follower spacecraft. The proposed control law is simulated for the case of the leader and 
follower spacecraft pair and is shown to yield global, asymptotic convergence of the relative 
position tracking errors. 


I. Introduction 

The Lagrange points of the Sun-Earth system have been exploited as key locations for space-based 
astronomical observation stations. 1,2 These locations are equilibrium positions in the restricted three body 
problem (RTBP), see Figure 1 (a) for details. The first three Lagrange points in the RTBP (labeled as Li, L2, 
and L3) are points that are collinear with the two primary masses (Sun and Earth). The last two Lagrange 
points in the RTBP (labeled as L4 and L5) are equilibrium points such that each of these points combined 
with the two primary masses yields an equilateral triangle. Each of the five equilibrium positions can host 
a spacecraft for an indefinite time period. A benefit of using a Lagrange point observation station is that 
spacecraft near these points obtain nearly an unobstructed view of the galaxy. Furthermore, missions near 
the Lagrange points are sufficiently far from the Earth, such that environmental effects (e.g., atmospheric 
and geomagnetic forces) do not affect spacecraft dynamics. 

Future space missions, that intend to utilize the L2 Lagrange point as the location for deep-space 
observations and/or interstellar communication have the advantage that solar influences on the spacecraft 
are minimal and space observations can be conducted on a frequent basis. In contrast, spacecraft that are to 
perform the same types of missions in either Sun-synchronous or low Earth orbits about the Earth are not 
suitable because these orbits expose the spacecraft to harsh physical conditions (e.g., gravitational and/or 
atmospheric disturbances, space debris, etc.). 

Recently, the European Space Agency has proposed the Darwin space mission, 4 which is to be deployed 
near the L2 Lagrange point where it will search for life in the universe and investigate the evolution of 
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galaxies. Scheduled for launch in 2014, the Darwin space mission will utilize six spacecraft to cooperatively 
work together in order to search nearby planets for traces of life, in the form of infrared radiation. 

An emerging technology to enhance space-based imaging/interferometry missions is spacecraft formation 
flying (SFF). SFF enhances space mission performance by distributing mission tasks, which are usually 
conducted by a monolithic spacecraft, to many small spacecraft. Thus, future space missions near the Sun- 
Earth Lagrange points can greatly benefit from SFF. However, to effectively utilize this new technology for 
space missions near the Sun-Earth Lagrange points requires proper design of spacecraft formations and for 
each spacecraft in the formation to be precisely controlled to maintain a meaningful baseline. 

Current spacecraft trajectory designs near the Sun-Earth Lagrange points consist of computing periodic 
trajectories in the form of Lyapunov and Halo orbits around a Lagrange point. 5-7 Unfortunately, these 
designs are used specifically to provide desired reference trajectories only for single spacecraft missions. In 
the current literature for SFF near a Lagrange point, a leader spacecraft is placed on a periodic orbit, e.g., a 
Halo orbit, around a Sun-Earth Lagrange point and a follower spacecraft is placed near this periodic orbit 
and a reference trajectory of the follower spacecraft relative to the leader spacecraft is designed. In Ref. 8, 
reference trajectories for follower spacecraft are computed using classical orbital elements, which result in 
bounded orbits around the leader spacecraft on a periodic orbit. In Ref. 9, feedback control is utilized to 
produce reference trajectories for follower spacecraft. In addition, Ref. 10 provides a method of generating 
reference trajectories for follower spacecraft using a numerical method, where the resulting trajectories are 
quasi-periodic. 

In this paper, we develop a leader-follower spacecraft formation, where the leader spacecraft is on a 
periodic, Halo orbit around the L 2 Lagrange point in the Sun-Earth system and the follower spacecraft is 
to track a desired relative trajectory. Specifically, we first develop the dynamics of the follower spacecraft 
relative to the leader spacecraft. Next, in the spirit of Ref. 10, we design a desired quasi-periodic relative 
trajectory for the follower spacecraft. In contrast to Ref. 10, our trajectory design exploits the analytical 
properties of the quasi-periodic relative trajectories to characterize spacecraft formations using a set of 
parameters. Finally, we develop an adaptive full-state feedback control algorithm to enable the follower 
spacecraft to track this desired quasi-periodic relative trajectory. 

This paper is organized as follows. Section II develops the mathematical model for the follower spacecraft 
relative to the leader spacecraft. Section III describes a method of generating follower spacecraft trajectories 
relative to the leader spacecraft orbit to create a spacecraft formation. Section IV formulates a trajectory 
tracking control problem. Section V uses a Lyapunov-based approach to design a full-state feedback con- 
trol law and a parameter update algorithm, which facilitate the tracking of given reference trajectories in 
the presence of unknown follower spacecraft mass. Illustrative simulations are included in Section VI to 
demonstrate the efficacy of the proposed trajectory generation and control design schemes. Finally, some 
concluding remarks are given in Section VII. 

II. System Model 

In this section, we develop a nonlinear model characterizing the position dynamics of the follower space- 
craft relative to the leader spacecraft near the L 2 Lagrange point in the Sun-Earth system. Referring to 
Figure 1, we assume that the Earth and the Sun rotate in a circular orbit around the Sun-Earth system 
barycenter (center of mass) with a constant angular speed to. In addition, we attach an inertial coordi- 
nate system {X,Y,Z} to the Sun-Earth system barycenter and a rotating, right-handed coordinate frame 
{xl 2 1 Ul 2 > z l 2 } t° the L 2 Lagrange point with the XL 2 -axis pointing along the direction from the Sun to the 
Earth, the ZL 2 -axis pointing along the orbital angular momentum of the Sun-Earth system, and the j/L 2 -axis 
being mutually perpendicular to the xl 2 and zl 2 axes, and pointing in the direction that completes the 
right-handed coordinate frame. 
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A. Dynamics of a Spacecraft Relative to the L2 Lagrange Point 

In order to describe the dynamics of a spacecraft formation near the L2 Lagrange point, we must first 
describe the dynamics of a spacecraft relative to the L 2 Lagrange point. To do so, let q(t) £) [x y z] T £ R 3 
denote the position vector from the spacecraft to the L2 Lagrange point, expressed in the {xl 2 , y\,. 1 . zl 2 } 
coordinate frame. In addition, let i?s->s(i) £ R 3 and R^^ s (t) £ R 3 denote the position vectors from the Sun 
and Earth, respectively, to the spacecraft. Finally, let R,j J2 , Re, and f?s denote the distances between the 
Sun-Earth system barycenter and the L2 Lagrange point, the Earth, and the Sun, respectively. Then, the 
mathematical model describing the position of a spacecraft relative to the L2 Lagrange point is given by 11 


mq + Cq + N(q , s) = u, 


(1) 


where m is the mass of the spacecraft, C is a Coriolis-like matrix defined as C = 2 mu> 
nonlinear term consisting of gravitational effects and inertial forces 
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and u(t) £ R 3 is the thrust control input to the spacecraft. Furthermore, the constants /ie and /ts in the 
definition of N are defined as he = GMe and /ig = GM$, respectively, where G is the universal gravitational 
constant, Me is the mass of the Earth, and Ms is the mass of the Sun. 


B. Halo Orbit Trajectory 

In this subsection, we describe a method to generate thrust-free, periodic trajectories around the L2 
Lagrange point in the form of Halo orbits. We present a succinct overview of a numerical algorithm to 
generate these periodic trajectories. Additional details on the generation of these periodic trajectories can 
be found in Refs. 5 - 7 . 

One numerical method' of generating thrust-free periodic orbits around the L2 Lagrange point in the 
Sun-Earth system involves finding a proper set of position and velocity initial conditions to propagate the 
spacecraft dynamics of ( 1 ), with the control thrust u set to zero. First, the Poincare-Lindstedt method is 
used to find a high order analytic approximation to a periodic trajectory in the neighborhood of the L2 
Lagrange point. Next, the initial conditions, based on the Poincare-Lindstedt method, are used as an initial 
seed in a numerical algorithm to find a better set of initial conditions leading to a periodic trajectory. This 
numerical algorithm applies a Taylor series expansion to the spacecraft states with respect to the initial 
conditions and time and truncates higher order terms, such that for Halo orbits the result is a set of 3 linear 
equations with 4 unknown variables. Families of orbits can be characterized by fixing one of the unknown 
variables so that the result gives an equal number of equations to unknowns. Solving the aforementioned 
linear matrix equation and using the result to update the previous set of initial conditions provide a new 
initial condition guess. 

The spacecraft dynamics are then propagated using the new updated set of initial conditions to verify 
trajectory periodicity. If the trajectory is sufficiently close to being periodic, then the initial conditions can 
be used for further simulation, else the above numerical algorithm is used to solve for a new set of initial 
conditions. Since the collinear Lagrange points are inherently unstable , 7 long-term propagation of spacecraft 
dynamics using the initial conditions obtained in the above manner is futile. However, by exploiting the 
symmetry property of Halo orbits (see below), we can artificially obtain a periodic orbit by computing 
trajectory information during half of a period and reusing this trajectory data throughout other simulations. 
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Halo orbits are classified as periodic trajectories that are symmetric with respect to the {xl 2 , Zl 2 } plane 
(i.e., j/l 2 = 0), and are not confined to be in the orbital plane of the Sun and Earth. Halo orbits have the 
distinguishing characteristic that their projections on the {?/l 2 ,~l 2 } plane are curves that resemble a Halo. 
In this paper, we let quit) = [#h(£) Unit) Zu{t)] T € R 3 denote the position vector from a point on a Halo 
orbit to the L 2 Lagrange point, expressed in the {xl 2 , Vl 2 , z l 2 } coordinate frame. An initial seed for the 
numerical algorithm of Ref. 7 consists of a spacecraft starting on the {xe 2 , <Zl 2 } plane with a nonzero initial 
Ul 2 and zl 2 velocity (i.e., <?h(0) = [xh(0) 0 Zh(0)] and gn(0) = [0 yn(0) id(0)] )• Updates to the 

initial Xl 2 position and y^ 2 velocity contribute to finding a closed periodic trajectory. In addition, the initial 
Zl 2 position determines the size of the Halo orbit. Figure 1(b) shows a typical Halo orbit trajectory around 
the L 2 Lagrange point. 

In this paper, we use Halo orbits as the reference trajectory for the leader spacecraft. The control design 
framework of Ref. 11 can be employed to ensure that the spacecraft dynamics of (1) tracks a Halo orbit 
reference trajectory. In a subsequent subsection, we will describe the dynamics of the follower spacecraft 
relative to the leader spacecraft on the Halo orbit. Finally, we denote Rs^u(t) £ R 3 and J?e->h(£) G R 3 as 
the position vectors from the Sun and the Earth, respectively, to the Halo orbit. 

Remark II. 1 The Halo orbit trajectory satisfies the spacecraft dynamics of (1) under the condition that the 
spacecraft control input is zero. Moreover, we express the leader spacecraft dynamics on the Halo orbit as 

mm + Cqu + iV(g H , H) = 0. (2) 

We note that the Halo orbit is a periodic trajectory with a frequency denoted as lo h- 

C. Follower Spacecraft Dynamics 

In this subsection, we describe the dynamics of the follower spacecraft relative to the leader spacecraft 
tracking a no-thrust, periodic Halo orbit trajectory gn without deviating from this orbit for all time. To 
describe the dynamics of the follower spacecraft, we express the position vector of the follower spacecraft 
relative to the L 2 Lagrange point in the coordinate frame {.tl 2 , 2/l 2 , -l 2 } as gf L (t) = [xf L yf Lo 2 f L ] £ 

R 3 . In addition, we denote Rs^ S! (t) £ R 3 and RE^ Sf (t) £ R 3 as the position vectors from the Sun and 
Earth, respectively, to the follower spacecraft. Using (1), the follower spacecraft dynamics relative to the L 2 
Lagrange point can be expressed as 

Wf« L2 + Cf5f L2 + Af L2 (q fh2 , Sf ) = u f , (3) 

A 0-10 

where mt is the mass of the follower spacecraft, Cf is a Coriolis-like matrix defined as C{ = 2m,fU> 1 0 0 

0 0 0 

IVf L is a nonlinear term consisting of gravitational effects and inertial forces defined as Nf h2 / t I ^-N(q{ L2 , Sf), 
and Uf(t) £ R 3 is the thrust control input to the follower spacecraft. 

Next, we define the relative position between the follower and the leader spacecraft qf(t) £ R 3 as qt = qi Lo — 
<7h- To obtain the dynamics of the follower spacecraft relative to the leader spacecraft, we differentiate qf 
with respect to time twice and multiply both sides of the resulting equation by mt to produce 

m f q f = m{q flj2 - rn f g H . (4) 

Next, we solve for gn hr (2), multiply the resulting equation by rrif, and substitute the result into (4) to yield 

m f g'f + C f q { + IVf(gf) = u { , (5) 

where (3) has been used. Note that Nf is a nonlinear term defined as Nf = Nf L (gf L , Sf ) — A r n(gH,H), where 
Ah is defined as N H ^^A(g H , H). 
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Remark II. 2 The Coriolis matrix C{ satisfies the skew symmetric property of x T CfX = 0, Va; £ R 3 . 

Remark II. 3 The left-hand side of (5) produces an affine parameterization mtqt + Cfqf + Nf(qt) = 
Y(qf, gf , qtfirif, where m{ is the unknown, constant mass of the follower spacecraft and Y(-) £ l 3 is a 
regression matrix defined as 
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III. Spacecraft Formation Design 

In this section, we exploit Ref. 10 to develop a method of designing reference trajectories for the follower 
spacecraft relative to the leader spacecraft on the Halo orbit trajectory. Specifically, we present a method of 
designing quasi-periodic orbits around a nominal Halo orbit. These quasi-periodic orbits will be used as the 
desired trajectories for the follower spacecraft. Furthermore, we will exploit special characteristics of these 
quasi-periodic orbits to parameterize spacecraft formations about the leader spacecraft on the Halo orbit. 

We begin by expressing the relative position dynamics of (5) in a state-space form, i.e., let x\ (t) £ R 3 be 
defined as x\ =qt and Xi (t) £ R 3 be defined as Xi =qf. Then (5) can be written as 


£2 

-mf 1 (C f x 2 + iVf(xi)) 

where Xf (f) ^ \xf xJ] T £ R 6 and we assume that Uf = 0, Vf > 0. Next, we linearize the nonlinear terms 
on the right hand side of (7), in the neighborhood of Xf = 0, to obtain 

Xf = AXf, (8) 

where A(t) £ R 6x6 is a time varying matrix with elements that are periodic with time. It is defined as 

A _i d-w/cxpi I fi„ , where 0s is the 3x3 zero matrix, Is is the 3x3 identity matrix, and 

m f dX! pi=o ”l J 

dN dxi^ 1^1=0 3S 3x3 Jacobian matrix of iVf(a;i) evaluated at x\ = 0. In addition, the period of oscillation 
of A is the same as the period of the nominal Halo orbit, i.e., A is periodic with a frequency wh. Note that the 
time dependence of A characterizes the dynamics resulting from the linearization of (7) as a nonautonomous, 
linear differential equation with a periodic A matrix. Consequently, we employ Floquet theory 12 to transform 
(8) into an autonomous, linear differential equation so as to facilitate an explicit solution of (8). 

We begin by introducing the notion of a fundamental matrix 12 of (8) denoted as ip(t) £ R 6x6 . Next, we 
denote the Halo orbit period as Th- Using Floquet theory, we utilize the transformation 

W = PYf, Yf = p-'Xf, (9) 

where Yf(t) £ M 6 is a vector composed of the transformed state Xf and P(t) £ R 6x6 is a matrix with elements 
that are periodic with time, 12 to transform the nonautonomous differential equation of (8) into 

Y { = BYf, (10) 

where B £ R 6x6 is a constant matrix. Following Ref. 12, the B matrix can be computed using tp and Th as 
follows B = log (<y2 _1 (0)(p(TH)), where the log function denotes the logarithm of a matrix. Furthermore, 
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the P matrix can be computed using p and B as follows P(t) = ip(t)e~ Bt . Note that the P matrix is 
nonsingular Vt G R, such that the transformation of (9) is unique. 13 

The autonomous, linear differential equation of (10) is equivalent to (8) in the transformed set of coor- 
dinates. Furthermore, the eigenvalues of the B matrix are denoted as the characteristic exponents, 13 which 
describe the stability characteristics of any trajectory that is sufficiently near the nominal Halo orbit. It 
is observed in Ref. 9 that direct computation of the eigenvalues of B results in a pair of hyperbolic eigen- 
values, a pair of zero eigenvalues, and a pair of nonzero, pure, imaginary eigenvalues. We denote the pair 
of hyperbolic eigenvalues as Ah! and Ah 2 and the frequency corresponding to the nonzero, pure, imaginary 
eigenvalues as ujq . Next, we perform a coordinate transformation of the form 


Y f = TZ { , 


( 11 ) 


where Zf(t) g K 6 is a vector composed of the transformed state Yp and T € R 6x6 is a time independent, 
linear transformation matrix, which transforms the B matrix into a modal matrix form, denoted by O G R 6x6 

0 1 1 I" 0 1 1 [" 0 1 

0 0 — AhjAp^ (Ahi + Ah 2 ) ~ lo ‘q ^ 

= nz { . 

Now, it is trivial to obtain the following solution of Zf analytically 


defined as fl^diag 


• . Then (10) is transformed into 


Y f = 


Yfi(0) + Z{ 2 (0)t Z{ 2 (0) 


-A h2 Zf 3 (0) + Yf 4 (0) ~ Xh t + A hl Zf 3 (0) - Z f4 (0) ^ t 

Ahi Ah 2 Ahi Ah 2 


— Ah 2 Yf 3 (0) + Yf 4 (0) 

Ahi ^h 2 


A hl e Ahl< + 


A hl Y f3 (0)-Z f4 (0) 

Ahi ^h 2 


A ha e Ah2 * D cos(wQt + </>) - Dloq sin(wQt + </>) 


,( 12 ) 


where Zf i (0), i = 1, . . . , 6, denotes the i th initial condition of the vector Zf and D, <fi G R are parameters 
that characterize size, location, and shape of the relative trajectory around the nominal Halo orbit. Eq. 
(12) reveals that the general solution of Zf may not be periodic for arbitrary initial conditions. However, 
by properly choosing the initial condition Yf(0) the terms corresponding to the pair of zero eigenvalues and 
the hyperbolic eigenvalues that produce unstable and/or asymptotically stable motion can be eliminated, 
thus resulting in periodic motion for Zf. The remaining periodic terms in (12) allow the trajectory designer 
freedom to choose the parameters Zf, (0), D , and </> to satisfy mission specifications. 

To compute the follower spacecraft trajectory relative to the nominal Halo orbit requires transformation 
from Zf — > Xf in the form of 

X t = PTZf, (13) 

where (9) and (11) have been used. Note that the P matrix is composed of elements which are periodic with 
respect to time, with frequency wh, whereas the solution to Zf is composed of elements which are periodic 
with respect to time, with frequency ujq. Consequently, the solution of Xf is a trajectory with two frequency 
components u > h and ujq. It is observed that these frequencies ujq and u > h are linearly independent, i.e., the 
condition ai^Q + 02 Wh = 0, a* G Z, i = 1, 2, where Z is the set of integers, holds only for a, = 0, i = 1, 2 (see 
Ref. 14 for details on linearly independent frequencies). Such a trajectory containing linearly independent 
frequency components are termed as quasi-periodic trajectories (see Ref. 14 for details on quasi-periodic 
functions). Thus, the Xf trajectory has the characteristic of being quasi-periodic. Finally, we utilize Xf as 
the desired trajectory of the follower spacecraft relative to the Halo orbit qd f {t) G R 3 , i.e., [<?J f = Xf. 

Remark III.l To facilitate subsequent illustrative examples, we approximate the Halo orbit and the P 
matrix using Fourier series approximations. Since both qu and P are periodic with the same period, the 
resulting Fourier series approximations are convergent to the actual forms of qu and P. To compute the 
time derivatives of qn and P, we analytically differentiate the Fourier series approximations with respect to 
time. Thus, it follows that qd { and its time derivatives, viz., qd { and qd { or equivalently Xf, are computed 
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using qn, P, and Zf, and their time derivatives, i.e., 


X t = PTZf + PTZf, 


(14) 


where (13) has been used. 

TV. Trajectory Tracking Problem Formulation 

In this section, we formulate a control design problem such that the follower spacecraft relative position 
qi tracks a desired relative position trajectory q^, i.e., lim qf(t) — q<\ s (t) = 0. The effectiveness of this control 

t— >oo 

objective is quantified through the definition of a position tracking error eft) £ R 3 as 

e = ®-< 7 d f - (15) 

The goal is to construct a control algorithm that obtains the aforementioned tracking result in the 
presence of the unknown constant follower spacecraft mass mf. We assume that the position and velocity 
measurements (i.e., qt and qi) of the follower spacecraft relative to the leader spacecraft on a nominal Halo 
orbit are available for feedback. 

To facilitate the control development, we assume that the desired trajectory q^ { and its first two time 
derivatives are bounded functions of time. Next, we define the follower spacecraft mass estimation error 
rhf(t) £ R as 

rhf £( rhf — mf, (16) 

where rhf ft) € R is the follower spacecraft mass estimate. 

V. Adaptive Position Tracking Controller 

In this section, we design an adaptive feedback control law that asymptotically tracks a pre-specified 
follower spacecraft relative position trajectory, in the presence of the unknown constant follower spacecraft 
mass mf. In order to state the main result of this section, we define the following notation. A filter tracking 
error variable rft) £ R 3 is defined as 

r £) e + ae, (17) 

where a £ R 3x3 is a constant, diagonal, positive-definite, control gain matrix. In addition, an augmented 
error variable is defined as qft) = \r T e T ] £ R 6 and a positive constant A is defined as A = min |A m i n {A'} , 

Amin {KpOi}}, where A m i n { - } denotes the minimum eigenvalue of a matrix and K,K p £ R 3x3 are constant, 
diagonal, positive-definite matrices. Next, we solve for e in (17) to produce 

e = r — ae. (18) 

Finally, we define a new regression matrix Id(-) £ R 3 as Vd(’) = T (£i , £ 2 , qf), where the linear parameteri- 
zation of Remark II. 3 has been used with £1 = q^ f — ae and £2 = q d f ~ ae, in the definition of (6). 

Theorem V.l Let K,K p £ R 3x3 be constant, diagonal, positive- definite matrices and T £ R be a positive 
constant. Then, the adaptive control law 


U{ = TdWf — K p e — Kr, rhf = —TYffr, 


(19) 


ensures global asymptotic convergence of the position and velocity tracking errors as delineated by lime(t), 

t — >-oo 

e(f) = 0. 
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Proof. We begin by rewriting the follower spacecraft relative position dynamics (5) in terms of the 
filtered tracking error variable (17). To this end, differentiating (17) with respect to time, multiplying both 
sides of the resulting equation by irif, using e = i:jf — q df from (15), substituting for mfijf from (5), and 
rearranging terms yield 

mff = — TOf(<jd f - cue) - Cf(j{ - N{(qf) + Uf. (20) 

Next, we expand (17) by noting that e = qt — q d[ . Then solving for qt, substituting the result into (20), and 
rearranging terms, we get 

mfr = — (g'd f — ae ) — Cf (gd f — cue) — Nf(qf) — Cfr + Uf = — Y d mt — Cfr + Uf, (21) 

where the definition of Yd has been used. Eq. (21) characterizes the open-loop dynamics of r. Now, 
substituting Uf of (19) into (21) results in the following closed-loop dynamics for r 

mff = Y d m { — K p e — Kr — Cfr, (22) 

where the definition of (16) has been used. Finally, note that differentiating (16) with respect to time and 
using rhf of (19), produce the closed-loop dynamics for the spacecraft mass estimation error 

fhf = -TY^r. (23) 

Now, we utilize the error systems of (22) and (23) along with the positive-definite, candidate Lyapunov 
function defined by V = \mfr T r + ^ e T K p e + ^ttt. 2 , to prove the above stability result for the position and 
velocity tracking errors. Specifically, differentiating V with respect to time and substituting the closed-loop 
dynamics of (18) and ( 22 ) into the result, we obtain 

V = —r T Kr— e T K p ae < — A|| 77 H 2 < 0, (24) 

where the property of Remark II. 2, (23), and the definitions of 77 and A have been used. 

Since V is a non-negative function and V is a negative semi-definite function, V is a non-increasing 
function. Thus V(t) € £00 as described by V(r(t),e(t),fhf(t)) < V(r(0), e(0), fhf (0)), t > 0. Using standard 
signal chasing arguments, all signals in the closed-loop system can now be shown to be bounded. Using (18) 
and ( 22 ) along with the boundedness of all signals in the closed-loop system, we now conclude that 77 G Coq. 

POO 

Solving the differential inequality of (24) results in U(0) — V(oo) > A / || 77 (t)|| 2 c?t. 

J 0 

Since V (t) is bounded, t > 0, we conclude that 77 (f) € £00 D ^ 2 , f > 0. Finally, using Barbalat’s Lemma , 15 
we conclude that lim 77 (f) = 0. Using the definitions of r and 77 , lim 77 (f) = 0, and Lemma 1.6 of Ref. 15, 

t—* OO t—> OO 

yield the result of Theorem V.l. □ 


VI. Simulation Results 

In this section, we present illustrative examples that incorporate the algorithms presented in Sections 
III and V. Specifically, we provide details on computing the quasi-periodic trajectories described in Section 
III. Next, we provide a simulation of the follower spacecraft relative dynamics (5), utilizing the control and 
adaptation laws of (19) so that the follower spacecraft tracks a desired quasi-periodic trajectory relative to 
a nominal Halo orbit. 

In all simulations, we employ the Sun-Earth system circular orbit parameters: 7, 16 G = 6.671 x 10~ n , 

u = 2.73774795629 x 10" 3 ^, M s = 1.9891 x 10 3 °kg, M E = 5.974 xl0 24 kg, 1 AU = 1.496 x 10 8 km, and 
i?L 2 = 1.010033599267463 AtJ, where 1 AU stands for 1 Astronomical Unit denoting the distance between 
the Sun and the Earth. Furthermore, we consider that the follower spacecraft has a mass of rrif = 1000kg. 
Finally, the distances f?s and Re can be computed as i?s = m e -m s x anc ^ = m e -m s x 
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A. Quasi-Periodic Trajectory Generation 

Applying the numerical algorithm presented in Subsection B results in a family of initial conditions for 
the Halo orbit from which we have selected the following initial condition gn(0) = [—2.61921376240742 0 
-0.13648677396294] x 10 5 km and q H ( 0) = [0 4.21353617291110 0] xl0 3 |^. In addition, the Halo orbit 

period is determined to be Th = 1.135225027876099 x 10 3 day. Figure 1(b) shows the Halo orbit relative to 
the L 2 Lagrange point and its projections onto the {a;L 2 ,yL 2 }, { x l 2 , z l 2 }, and {2 /l 2 >~l 2 } planes. In addition, 
we utilized 25 terms of a Fourier series to approximate the Halo orbit trajectory q h. The fundamental matrix 
tp described in Section III is numerically computed using A(t) as follows p = A{t)p, y>(0) = Ie, Vi € [0,Th]. 
Thus, using P(t) = p(t)e~ m , P is numerically computed Vi £ [0,Th]. Next, we compute a Fourier series 
approximation of P, where we retain 25 terms of the series approximation. This, along with the analytic 
expression for Z{ is used to compute and its time derivatives analytically from (14). 

To show the resulting trajectories of q a f , given different numerical values for parameters Zf 1 (0), D, and 
0, we simulated using a parameter set: Zf x ( 0) = 0, D = 0.0001, and 0 = 0 rad. By computing the 
eigenvalues of the B matrix, we determined ujq = 6.286301816644046 x 10 -5 ^d_. Figure 2(a) shows the 

quasi-periodic trajectory relative to the nominal Halo orbit for parameter values of 0 = 0, 0 = and 
0 =t^. Figure 2(a) illustrates that changes in 0 denote changes in the initial position of the spacecraft along 
a given quasi-periodic trajectory. Next, we simulated q^ f using a parameter set: Z{ 1 ( 0) = 0, D = 0.0002, 
and 0 = 0 rad. Figure 2(b) shows the desired quasi-periodic trajectory relative to the nominal Halo orbit. 
Note that the parameter D determines the size and shape of the desired quasi-periodic trajectory relative 
to the nominal Halo orbit. We also simulated gd f using a parameter set: Z{ 1 ( 0) = 0.0001, D = 0, and 0 = 0 
rad. For this parameter set, Figure 2(c) shows a periodic trajectory relative to the nominal Halo orbit with 
the same period as wh. Finally, we simulated q<i { using a parameter set: Zf 1 ( 0) = 0.0001, D = 0.0001, and 
0 = 0 rad. For this parameter set, Figure 2(d) shows the quasi-periodic trajectory relative to the nominal 
Halo orbit. 

B. Adaptive Full-State Feedback Control of Follower Spacecraft 

The adaptive control law of (19) was simulated for the follower spacecraft dynamics relative to the leader 
spacecraft on a nominal Halo orbit (5). When tracking desired quasi-periodic trajectories, we initialized the 
follower spacecraft with the set of initial conditions given as <7f(0) = [—2.61921376240742 — 2.57780484325713 
-0.13648677396294] x 10 5 km and ®(0) = [-0.01469110370264 4.21353617291110 - 0.01469092330256] x 
10 3 ^|. The control and adaptation gains are obtained through trial and error in order to obtain good 
performance for the tracking error response. The following resulting gains were used in this simulation 
K = diag (1,1,1) x 1.499 x 10, K p = diag (1,1,1) x 5.475 x 10 3 , a = diag (1,1,1) x 8.213 x 10~ 2 , and 
T = 8.888 x 10 4 . In addition, the follower spacecraft mass parameter estimate was initialized to mf(0) = 600 
kg. A simulation of the follower spacecraft tracking the desired quasi-periodic trajectory of Figure 2(a) is 
performed. The trajectory qt is shown in Figures 2(e) and 2(f). Figure 3 shows the position tracking error 
e and the velocity tracking error e. The control input Uf is shown in Figure 4(a). Finally, the follower 
spacecraft mass estimate rh{ is shown in Figure 4(b). 

VII. Conclusion 

In this paper, we designed desired quasi-periodic trajectories for the follower spacecraft relative to the 
leader spacecraft on the Halo orbit. The size, location, and shape of these trajectories were characterized by 
a set of parameters. Illustrative simulations were performed to show these parameter characteristics. Next, 
a Lyapunov design was used to develop an adaptive full-state feedback controller, which yielded global, 
asymptotic convergence of the relative position tracking errors. Simulation results were presented to show 
good trajectory tracking. 
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Figure 1. (a) Sun-Earth system schematic diagram and (b) Halo orbit trajectory of the leader spacecraft 

relative to the L 2 Lagrange point 
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Figure 2. Trajectory of the follower spacecraft relative to the nominal Halo orbit using: (a) Zf x (0) = 0, D = 
0.0001,0 = 0,0 = f , and 0 = f , (b) Z fl (0) = 0, D = 0.0002,0 = 0, (c) Z fl (0) = 0.0001, D = 0, 0 = 0, (d) 
Zf x (0) = 0.0001, D = 0.0001, 0 = 0, (e) normal view, and (f) zoomed in view 
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Figure 3. (a) Position tracking error and (b) velocity tracking error 
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Figure 4. Follower spacecraft (a) control input and (b) mass parameter estimate 
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